Nonlinear acoustic and microwave absorption in disordered semiconductors 
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Nonlinear hopping absorption of ultrasound and electromagnetic waves in amorphous and doped 
semiconductors is considered. It is shown that even at low amplitudes of the electric (or acoustic) 
field the nonlinear corrections to the relaxational absorption appear anomalously large. The physical 
reason for such behavior is that the nonlinear contribution is dominated by a small group of close 
impurity pairs having one electron per pair. Since the group is small, it is strongly influenced by 
the fleld. An external magnetic field strongly influences the absorption by changing the overlap 
between the pair components' wave functions. It is important that the influence is substantially 
different for the linear and nonlinear contributions. This property provides an additional tool to 
extract nonlinear effects. 

PACS numbers; 72.20.E, 71.23.A 
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I. INTRODUCTION 

The subject of this article is nonlinear microwave and 
acoustic properties of amorphous semiconductors and 
lightly doped crystalline semiconductors in the regime of 
hopping conductance. We are interested in absorption to 
electron transitions between localized states associated 
with defects or impurity atoms. We consider the case 
where the absorption is due to electron hopping within 
pairs of neighboring defects containing one electron per 
pair. The distance between the centers within the pair 
must be small enough to allow tunneling, while the dis- 
tance to other impurities should be large enough to pre- 
vent tunneling to impurities outside the pair. In a weakly 
doped semiconductor we can expect these pairs to be rel- 
atively rare, and triplets of the same kind will be even 
less likely. Thus a natural basis to treat the problem is 
the so-called two level approximation according to which 
only the lowest energy level of each of two neighboring 
impurities are taken in consideration. This approach, as 
well as its range of applicability, was first discussed in 
detail by Pollack and Geballcld 

For brevity, in the following we shall discuss the case 
of acoustic attenuation and then specify what changes in 
the formulae should be introduced to allow for electro- 
magnetic absorption. 

An external AC electric or acoustic field causes transi- 
tions between the electron states. Direct inter-level tran- 
sitions leading to absorption of quanta give rise to the 
so-called resonant absorption. For low intensities, the 
resonant contribution to the absosption coefficient of an 
acoustic wave can be expressed asa 
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where ai is a dimensionless coupling parameter, weakly 
dependent on temperature and frequency, s is the sound 
velocity, uj is the sound frequency, and T is the tempera- 
ture. We can always assume the relation hui/2kBT <C 1. 



ai will be specified later, see Eq. (p^). Defining hq 
as the Fermi function no{E) — [exp(£^/fcsT) + 1]^^, 
we can write the factor tauh(huj/2kBT) as no(?ia;/2) — 
no{~hLu/2). This can be recognized as the difference in 
the equilibrium population of the two levels of the pair 
with an energy splitting of huj/2 — {—fiuj/2) = hoj, that 
is, the pairs which can directly absorb a phonon. 

The relaxation absorption is due to a modulation of 
the electron inter-level spacing 2e by the AC field. Such 
a modulation leads to a periodic change of the occupation 
numbers of the two levels which lags in phase the varia- 
tion of e. This lag leads to the energy dissipation. In the 
linear regime the coefficient of relaxation absorption has 
been calculated asQ 
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where 012 ~ ai and tq represents a minimal relaxation 
time for e « fcsT. The physical meaning of tq, as well as 
estimates of this quantity, will be will be discussed later. 

Comparing Eq. (H) with Eq. (|l|) we conclude that the 
relaxation absorption always predominates at cjtq ^ 1. 
If WTO » 1 the ratio r(''<='')/r(''<=') « LOTQtimh{huj /2kBT) 
can be either greater or less than one under experimen- 
tally accessible conditions. 

For higher intensities this comparison is no longer 
valid, as both the resonant and the relaxation absorp- 
tion show strongly nonlinear behavior. For the resonant 
absorption, the nonlinearity is due to an equalization of 
the population numbers of the two electron states. From 
Eq. (|l|) we see that this leads to a strong reduction of 
the resonant absorption. Usually this suppression takes 
place at very low intensities.H 

The nonlinearity of the relaxation absorption is due 
to the following. In the limit of high intensities, when 
the perturbing potential is amplitude ^ fc^T, there will 
be times when 2e{t) < kBT and both states are almost 
equally occupied. During such part of the wave period no 
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transitions and thereby no absorption will occur. Con- 
sequently, the absorption coefhcient decreases with the 
sound amplitude. 

The linear relaxation absorption in semiconductors has 
previously been studied both theoreticallyul and experi- 
mentally (for a review, see Ref. ||), while-for the strongly 
nonlinear regime only theoretical resultaa were obtained. 
The reason why this type of nonlinearity has not been ob- 
served experimentally is probably some masking by other 
mechanisms leading to a nonlinear behavior. One of such 
mechanisms could be wave-induced ionization of injmurity 
atoms into the conduction band obser|Ved in InSbfl Con- 
sequently, to observe the mechanismo of nonlinear ab- 
sorption one should carefully choose the material. That 
does not seem to be an easy task, and we are not aware 
of any experiments of this type. 

In this article we will address the case of low intensi- 
ties when nonlinear effect manifest themselves as small 
corrections to the linear absorption. The main message 
is that already the lowest order corrections offer inter- 
esting information and anomalous effects worth study- 
ing. Hopefully, these effects will be pronounced within an 
experimentally achievable parameter space where other 
mechanisms of nonlinear behavior are still not important. 

The two-site approximation for semiconductors allows 
us to describe our system in terms of the so-called Two 
Level System (TLS) model. This moclel was first pro- 
posed independently by Anderson et alE and PhilipsB to 
explain the low temperature specific heat in amorphous 
materials, and has been successful in describing several 
other phenomena. The model has given a theoretical ex- 
planation for a surprising universality in the behavior of 
very different disordered materials at low temperatures. 

The authorsEI have previously used the TLS model to 
analyze nonlinear corrections to the absorption in both 
dielectric and metallic glasses. However, there are several 
important differences between the situations in glasses 
and in disordered or doped semiconductors. The main 
difference is due to the electric charge of the particle in- 
volved in the transitions. This introduces two modifica- 
tions to the earlier results. 

First, the electric charge gives the pair a dipole mo- 
ment proportional to the distance between the pairs, the 
coupling being proportional to the dipole moment's com- 
ponent along the direction of the electric field. This leads 
to a specific orientational dependence of the absorption 
of one pair, and after integration over all pairs it signifi- 
cantly infiuences the absorption. Second, we must expect 
our system to behave differently with the application of 
a magnetic field. More specifically, as the magnetic field 
leads to a stronger localization of the electrons, we must 
expect the absorption to decrease with an applied field. 
The effect of a magnetic field will therefore also be ana- 
lyzed, in the special cases of a weak or a strong field. The 
connection between the orientation of the dipole moment 
and the magnetic field direction also leads to a depen- 
dence of the absorption coefficient on the relative direc- 
tions of the radiation wave vector and the magnetic field. 



We will solve the problem for a magnetic field parallel or 
perpendicular to the wave vector of the radiation. 

The magnetic field dependence gives us the possibility 
to separate the relaxational absorption from other con- 
tributions. It will also be shown that the effect of a mag- 
netic field depends on the specific relaxation mechanism, 
thus providing us with a tool to further understand the 
relaxation processes for localized carriers. 

The paper is organized as follows. First, we give a short 
introduction to the theory based on the TLS model and 
how it can be used to solve the problem of relaxation 
absorption in glasses. Then we will fit the two-site ap- 
proximation for a semiconductor to the TLS model and 
show what modifications are needed for this. Finally, we 
will show how the effects of a magnetic field can readily 
be included in the analysis, and what effects to expect 
from this. 



II. TWO LEVEL SYSTEMS AND RELAXATION 
LOSSES IN GLASSES 

The TLS model deals with a particle (in our case, an 
electron) moving in a slightly asymmetric double well po- 
tential. It is assumed that only the two ground levels are 
accessible. The ground levels in the isolated wells are 
assumed to have a slight separation in energy, 2 A, and 
to be coupled via a tunneling energy overlap integral A. 
The Hamiltonian of such a system is traditionally written 
as 



TYo = A CT^ - A CT;i 



(3) 



where tXi are the Pauli matrices. 

Let us now apply an external periodic perturbing po- 
tential and study the power absorbed by a single pair, 
p(A,A). As is will be shown, this power is in general 
a non- monotonous function of its parameters, and there 
exists an "optimal" region which dominates absorption. 
What is important is that those regions are different for 
the linear absorption and for the nonlinear correction. 
Let us define S as the typical value of the inter-level 
splitting A which is important for the onset of nonlin- 
ear behavior. Estimates for the quantity S will be given 
in the discussion section. 

At huj <IC S one can employ the adiabatic approxi- 
mation and neglect time derivatives of the externalp&eld 
while solving the Schrodinger equation for the TLSJlj In 
t^ approximation we write the interaction Hamiltonian 



'Hj = Gzd cos ojt 



(4) 



ignoring possible off-diagonal items The quantity d is 
just the coupling constant between the field and the TLS. 
In the physics of low-temperature properties of glasses d 
is assumed to be a random quantity, uncorrelated with 
A and A. This assumption is generally not valid for the 
case of semiconductors. 
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In the case of a sound wave d = "fikuf^ , where is 

the deformational potential of the TLS and is the 
amphtude value of the deformation tensor. 

For semiconductors the value of A depends on the spa- 
tial separation of the wells. In the case of electromag- 
netic waves d = r]£o, where ij is the dipole moment of like. 
TLS while Sq is the amplitude of the electric fieldE'EJ 
For charged particles the dipole moment is proportional 
to the distance between the wells, and is thereby strongly 
correlated with the value of A. 

The total Hamiltonian of the TLS may now be written 

ai 



H = {A + d cos ujt) (jz — ^(Tx 



(5) 



The difference between the eigenvalues of this new Hamil- 
tonian, 2e(i), where 



e{t) = V(A + dcoscjt)2 A2. (6) 

is just the energy splitting of the TLS.I1I To characterize 
a TLS we need not only the energy spacing, but also the 
occupation numbers of the upper (n) and lower (1 — n) 
levels. The non-equilibrium occupation numbers can be 
found from the balance equationll3 



dn n — noit) 

H " T{t) 



(7) 



where no{t) is the adiabatic equilibrium occupation num- 
ber. It depends on the energy spacing e(t), temperature 
T and time t as 



noit) 



,2eit)/kBT 
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(8) 



The relaxation time T{t) is a function of the energy split- 
ting, the tunneling barrier and temperature, and also de- 
pends on the exact relaxation mechanism. 

The power abscpbed by a single TLS can be determined 
by the expressiont3 



p(A,A) = - / dtn{t)- e= — 
vy Jq dt uj 



(9) 



The contributions of individual TLS must be added 
and such a summation can be performed in a conven- 
tional way using the distribution function A^(A, A) of the 
random parameters A and A and replacing the deforma- 
tional potential 7^^ by its average value. 

To analyze the nonlinear absorption we use the exact 
periodic in time solution of Eq. (0) to-obtain the following 
result for the total absorbed power ,113 
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Unfortunately, this integral cannot be calculated an- 
alytically in the general case. Earlier discussions of 
the strongly nonlinear regime have considered situations 
where only very limited time or energy ranges contribute 
to the absorption. This corresponds to the asymptotic 
high intensity limits of the integral. Our approach has 
rather been expanding the integral in powers of the am- 
plitude of the modulation of the energy splitting, d, to 
get insights into the first onset of nonlinear effects. So 
we will concentrate on the regime of weak nonlinearity. 



III. SEMICONDUCTORS IN THE TLS 
FORMALISM 

We will now show how hopping in amorphous or doped 
crystalline semiconductors can be described by the TLS 
model. 

A TLS is naturally formed by a pair of nearest impu- 
rity centers having one electron per pair. To make the 
calculations for an individual pair we have to specify the 
coupling constant d and the relaxation time r. The latter 
is essentially dependent on the dominant mechanism of 
pair population relaxation. To sum over all the pairs we 
have to specify the proper distribution function for the 
parameters of the pairs containing one electron. 

We will concentrate on the case when the interaction 
has a dipole character and can be expressed in the form 



d{r,t) ^ do (r) cos ujt , do (r) = e£or{v ■ n) 



(11) 



where Eq is the amplitude of the effective electric field 
acting upon the electrons. It is just the amplitude of the 
local electric field created either by the external electric 
field, pur due to piezoelectric interaction with an acoustic 
wave.liJ v> — £/£ is the field polarization vector, r is the 
distance between the components of the pair, n = r/r is 
the pair direction vector. 

The relaxation time is strongly dependent on the par- 
ticular mechanism of interaction between localized elec- 
trons and phonons, see for example review Ref. ^. Here 
we will only quote the most important results. 



A. Relaxation time 

In general, the relaxation time can be expressed aJ 

2 
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A(r) 



2e \ coth{e/kBT) 



(12) 



where the exponent n, the factor and the function 
$n(a:) are dependent on the particular mechanism of in- 
teraction between localized electrons and phonons. The 
meaning of the energies Er and Ea will be made clear in 
a moment. 
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The energy dependence of can be easily under- 
stood. The phonon emitted or absorbed during a tran- 
sition between the electron levels must obviously have 
an energy equal to the energy splitting of the pair, 2e. 
The power by which it occurs in the formula ( p^ ) is de- 
termined by the product of the phonon density of states 
with the frequency dependence of the squared interac- 
tion matrix element. The factor coth(e/fcsr) is equal to 
2N^ + 1 where N,^ is the Planck function for Uui = 2e. 
Apart from a proportionality factor, this is the proba- 
bility of phonon emission, cx N^^ , plus the probability of 
absorption, cx {N^ + 1). The factor [A(r)/e]^ is a di- 
mensionless measure of the tunneling coupling between 
the bare states in the two wells. It can be seen that r 
has a minimum with respect to A when A = e, which is 
equivalent to A = 0. This condition defines the minimal 
relaxation time tq referred to in Eq. (|^). The minimal r 
corresponds to the symmetrical configuration when the 
bare energy levels at both sites are equal. To allow for 
time dependence of r one should substitute e = e{t) from 
Eq. (I). 

The expression (^2|) contains two specific energy scales, 
Ea = 2Tis/a, and Er = hs/r. The first scale is the en- 
ergy of a phonon having a wavelength of the order of the 
single-site localization length, a. Phonons with larger en- 
ergies produce rapidly oscillating fields which average out 
at the distance occupied by a localized electron. Conse- 
quently, strongly decays at e > Ea- The second 
scale corresponds to phonons with a wavelength of the 
order of the distance r between the components of the 
pair. If the deformation potentials of both components 
of the pair are the same, or if the main mechanism of 
the electron-phonon interaction is piezoelectric, then at 
2e <C Er both levels move synchronously, and no interac- 
tion occurs. The net interaction is hence proportional to 
some power of the ratio 2e/Er, see Rcf. || for a review. 

In most realistic cases one can assume e <^ Ea- How- 
ever, for T « 1 K the ratio x = 2e/Er can be either less 
or greater than 1, giving different types of behavior for 
^n{x)- We will concentrate on the case of x ^ 1 since 
this limiting case seems to be more easily accessible for 
experiments. In this regime, ^n{x) can be considered as 
constant. The validity of this approximation will be con- 
sidered in the Discussion. The quantities t„ (T) are listed 
in Rcf. I 

The most important feature of the relaxation for our 
problem is the energy dependence of the relaxation rate, 
namely, the power n. Under the above-mentioned con- 
ditions, n = 3 in the case of deformational interaction 
and n = 1 for piezoelectric interactions o As will be clear 
later, only the energies e <^ ksT are important for the 
anomalous nonlinear behavior, so one can approximate 
coth(e/fcBT) « ksT/e- In this way we arrive at the fol- 
lowing energy dependences of the relaxation rate: 



1 

— cx 



A^e° for the deformational interaction 
A^e^^ for the piezoelectric interaction 



(13) 



dent on e and thereby on time. This is the same as is the 
case in dielectric glasses. Apart from a constant and the 
dependency on the magnetic field, we can thus expect the 
same type of behavior from these two very different sys- 
tems. In the second case the e-dependence is the same as 
in mctaUic glasses, which has also been analyzed by the 
authorsO and has proven to be the source of a pronounced 
anomalous effect. In particular, in metallic glasses the 
lowest nonlinear contribution is proportional to the in- 
tensity to 3/2 rather than to the intensity squared as for 
dielectric glasses. 

In the absence of the magnetic field the energy over- 
lap integral A is related to the distance r between the 
sites of a pair simply as A = Aoexp"'"/", where Ao = 
(1 — 5) X me^/?i^K^ is of the order of the effective Bohr 
energy. Here m is the electron effective mass while n is 
the dielectric constant. A magnetic field will squeeze the 
electron wave function, and this effect will be strongest 
for the direction perpendicular to the field. This intro- 
duces an angular dependency to the localization length 
and thereby also to A. Following Ref. |l^, we will ana- 
lyze the limiting cases of weak (w) and strong (s) mag- 
netic field where the inffuence of magnetic field is weak 
or strong, respectively. The asymptotic expressions for 
the C(r) = -ln[A(r)/Ao] are the following: 



r r a sin 6 
a ^ 24A4 ' 
sin^ 9 |rcos( 



4A2 



an 



(14) 
(15) 



Here 6 is the angle between r and the direction of the 
magnetic field H, A = y/Tic/eH is the magnetic length, 
while an = h/y/2mEH is the characteristic localization 
length in the longitudinal direction, where Eh is the ion- 
ization energy of the ground state of the localized electron 
in the magnetic field. 



B. Pair distribution function 

The total absorption is given by a sum of the contribu- 
tions of the individual pairs. Hence we have to sum over 
the r, as well as over the individual energies of the elec- 
tron levels. The latter summation must take into account 
the correlation between the level occupation numbers due 
to Coulomb interaction. As shown in Rcf. 16 , the summa- 



tion over the energies can be split into integration over 
the pair center-of-gravity and over the bare inter-level 
spacing A. The first integration gives 2 A -I- e'^/nr since 
only the pairs with the center-of-gravity energy between 
the chemical potential /i 4- A and /i — A — e^/^cr have one 
electron per pair. As a result, the pair distribution fmic- 
tion can be expressed through the single-electron density 
of states g aal3. 



Note that in the first case the relaxation time is indcpcn- 



iV(A,r) 



g^V 
47r 



2A 



e 

Kr 



(16) 
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where V is the volume contributing to the absorption. 
This expression is vahd if for a typical hopping distance 
e^/nr ^ Ac where Ac in the width of the Coulomb 
gap in the single-electron density of states Il3 Inside the 
Coulomb gap, we would rather have to use the distribu- 
tion 

Calculations for both cases are similar. Note that the 
distribution function is isotropic. Under the conditions 
of interest to us the typical hopping distance r is small 
enough to let us neglect 2 A in comparison with e'^/nr 
in Eqs. (|l^) and (p^). Thus the distribution becomes 
A-independent, and we denote it as N{r). Using the no- 



tation presented above we can now write the |-cpupling 
constants ai, alpha2 ~ a from Eqs. |l|, and || asllZl 

where JC is the coupling constant of the piezoelectric in- 
teraction. The power of r^^ may vary for different type 
of interactions, depending on whether the interaction in- 
cludes the dipole moment of the pair. 

IV. CALCULATION OF ABSORPTION 

As a result of the previous considerations, the absorbed 
power can be expressed as 



P= dn 



dA 



drN{r) 



/.e 



-'0 



dtdt' e{r,t)e{r,t-t') ^xp 



QkgT cosh^ [e(r, t - t')/kBT] 



exp[ 



(19) 



Here the expression in the braces is just the power ab- 
sorbed by an individual pair, p(A, r). The energy split- 
ting 2e depends on r through the interaction potential 
(io(r) given by Eq. ( |ll| ) and through the tunneling split- 
ting A = Aq e~^^""^ Technically it is convenient to trans- 
form the integral from the set of variables r, n, to the 
variables A, 9, (j). Such a transform introduces the factor 



9A(r, n) 



dr 



Here the dimensionless function /(A, n) is given by the 
equation 



•'^ a d\nA a dC 



(20) 



where £ = ln(Ao/A), while rA(n) is the solution of the 
equation 



C(r-A,n) = £. 



In the simplest case of zero magnetic field / = 1, and 
we have a distribution very similar to those of glasses. 
The variable transform strongly simplifies the calcula- 
tions since the relaxation time is a function of the param- 
eters A and A and the quantity r\ is a weak (logarithmic) 
function of A. Consequently it can be extracted out of 
the integral over A, while replacing A in the expression 
for TA by its characteristic value. As a result, the inte- 
grations over A and A will remain the same as previously 
calculated for glasses, and only the angular integration 
and the dependence on the characteristic value of A are 
different. 



The following calculation procedure is similar to that 
of Ref. 17. The expression ([l9| ) will be expanded in pow- 
ers of the effective electric field, and the lowest correction 
will be compared with the linear result. Expanding the 
individual contributions in powers of do as 



p(A,r)=^pWdg(r) 



(22) 



k=2 



we notice the the coefficients p^'^^ depend only on the 
quantities A and A. Transforming the variables from r, n 
to A, n we can use the fact that r is a weak function of A 
and extract of the quantities proportional to the powers 
of r from the integral over A replacing 

r ~> rc{n) = rA(n) |a=a, 

where Ac is the characteristic value determined by the 
integrand over A. In a similar way, we replace /(A, n) — s- 
/c(n) = /(Ac,n). Finally we arrive at the expression 



(21) p = J2t=2 P^^^ with P^^^ ^ {aeSofhJk where 



^k=2 

f'OQ /"OC 

h = dA dAA-^p^''\A,A), 
Jo Jo 

Jk - / dn{vnff,N{r,){rJaf 



(23) 
(24) 



The quantities Ik are the same that enter the expressions 
for nonlinear absorption in glasses, and we quote them 
from Ref. ||. 

Parallel with Ik we can also estimate Ac, and thereby 
C. It can be shown that this value is only weakly depen- 
dent on the mechanism of absorption and on k. For all 
the mechanisms analyzed in this paper we have 
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£ = In ■ 



c^Ao 



(25) 



The value for cc is of the order one, and varies only with 
a factor less than two for the situations discussed in this 
paper. We will therefore treat C as independent of k for 
this discussion. 

In this paper we will present the results for the case 
of low frequencies, uiTn «C 1, for which the anomalous 
nonlinear behavior is most pronounced. In this case the 
linear in the intensity contribution, /q , is independent of 
r, and moreover, it is the same for any dependence T(e), 



/o (7rVl6) w w 0.62 w. 



(26) 



The case of deformational interaction, for which the 
parameter n in Eq. ( p^ is equal to 3, corresponds to the 
situation in dielectric glasses. In this casetl 



, /4 ~ 0.054 



(27) 



The piezoelectric interaction {n — 1) is similar to the 
case of metallic glasses. It can be shown that the energy 
dependence of the relaxation time leads to a divergence 
in the integration over A and A. To get a proper estimate 
one should cut off the integration at e ^ do. As a result, 
the leading nonlineac, contribution appears proportional 
to jdol'^ and equal toEl 



(28) 



We again expect C to be given by Eq. P5| . To calculate 
the angular integrals Jk we solve Eq. (|2l|) in a recursive 
way. For the case of weak magnetic fields, we obtain 



r^ = aC{l- C^a'^sm^ 9/2A\'^) 
fc= l + a'^C^cos^e/SX'^ . 



(29) 
(30) 



For the strong field limit, solving a quadratic equation, 
we obtain 



2A2 cos 9 

Tc — 7-5— 

an sin 9 
fc = (an /a) £ cos 



1 



1 



tan2 9C{aH/Xy - 1 



(31) 
(32) 



Here 9 is the angle between n and the direction of mag- 
netic field H. Hence, dn = d{cos9)d(f) where (p is the 
azimuthal angle between the projections of n and £0 
on the plane, perpendicular to H. After substitution 
of Eqs. (^-(32) into Eq. (|2^ ) the angular integrals are 
calculated directly. 



V. RESULTS 

To set the scale of nonlinear corrections let us start 
with the expression for the linear absorption in the ab- 
sence of magnetic field. Pa = Po(0). Both linear and 
nonlinear contributions are dependent on the relaxation 
mechanism of the relevant pairs. 



A. Deformational interaction between localized pairs 
and phonons 



We will first consider deformational interaction be- 
tween the localized pairs and thermal phonons. For the 
deformational mechanism, 



Po = (^V48) {Va^g^e^LuC^£^/K) 
In the absence of magnetic field we obtain 
P{0)-Po{0) 



P4(0) 



^o(O) 



0.26 



^wro(r) 



(33) 



(34) 



were we have introduced the dimensionless "field ampli- 
tude" 



F = eSoa/ksT . 



(35) 



In a weak magnetic field, the quadratic in magnetic field 
corrections arise both to the linear absorption and to the 
lowest nonlinear contribution. They can be expressed in 



a unified way as, cf. with Ref. 17 



ij^iH)^Po/m [i-cUa/xrc 



4 /.21 



(36) 



so the magnetic field produces corrections which are 
(X H^. The numerical factor depends on the direc- 
tion of the electric field £0 with respect to the magnetic 
field H. Its values are also different for the linear and 
nonlinear, contributions, Cw^ and Cw\ respectively. The 
values of Cw are shown in Table I. 



Direction 


((') 




^(3) 


H||5o 


0.2 
0.1 


« 0.29 
w 0.095 


« 0.25 
w 0.097 



Table I. Numerical coefficients entering the 
nonlinear contributions to the absorption. 

The decrease of attenuation in the magnetic field has 
the following physical reason. The presence of a mag- 
netic field squeezes the electron wave functions, and the 
overlap integrals between the components of the pair de- 
crease. Furthermore, the wave functions are squeezed 
mostly in the direction perpendicular to H. On the other 
hand, the coupling between the wave and the pair is max- 
imal if the pair dipole moment is parallel to £a. Thus the 
reduction of absorption is more pronounced for H _L n. 
For the nonlinear contribution the difference should be 
even stronger, as it includes higher orders of the dipole 
moment. 

In the limit of strong magnetic fields the results are 
even more interesting, as the functional dependency on 
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the magnetic field also varies with the different absorp- 
tion types. Yet we still see the same relative consider- 
ations as for the weak field limit. For a magnetic field 
parallel to the radiation polarization vector we get 



2„2 
H 



Pa Ca^ 



OC 



P4(0) 2 £a6 



(37) 



(38) 



For a perpendicular field the results show both a stronger 
dependency on H and on the order of the expansion in 
intensity. 



- 6^r^ln^# cx i/^^lnF, 



^'o(O) 



A2 
In— # 



(39) 



(40) 



B. Piezoelectric interaction between localized pairs 
and phonons 

We will now turn our attention to the piezoelectric in- 
teraction. This also interacts via a dipole moment, so 
the linear results are basically the same, apart from a 
coupling constant. However, there is a striking difference 
for the nonlinear contribution - similarly to the case of 
metallic glasses, B the integration over A and A in Eq. ([23 ) 
results in a term proportional to |(iop rather than dg. 
This is reflected in a change of the absorption depen- 
dencies both on the wave intensity and on the magnetic 
field. For this case, we restrict ourselves by order-of- 
magnitude estimates for numerical factors. Calculation 
of exact numbers would require a great amount of numer- 
ical work which would be inadequate to the accuracy of 
the initial model for the electron density of states. The 
results read as, 



P3(0) = [P-P(0)]/Po = c,„|F|, 



(41) 



where the estimate for Cm is 0.1. Note that in this case 
the expansion of the absorption in powers of intensity ap- 
pears non- analytical which implies relatively strong non- 
linearity. 

Since now the expansion in d of the nonlinear contribu- 
tion starts with |doP rather than from dg one can expect 
weaker magnetic field effects. This is indeed the case. 
The values of the numerical coefficients are shown in 
the Table I. 

For strong fields we get for parallel and perpendicular 
fields respectively 



V ^3(0) ^ 

pm J 



8A^ 
3 Ca^ ^ 

64 A^ 



(42) 



In:^ oci/-^/2lni/. (43) 



Let us recall that the above expressions are valid for the 
low frequency limit only, where ujTn{T) << 1. Similar 
calculations are possible for the high frequency limit, as 
well as for different pair distribution functions. In par- 
ticular, for the case of pronounced Coulomb gap the es- 
sential differences occur only in powers of r and thereby 
of C. Thus the influence of the magnetic field will be dif- 
ferent. The relation between linear and nonlinear results 
remain similar, apart from numerical factors. 



VI. DISCUSSION 



Let us discuss the relevance of the obtained results for 
realistic materials and situations. In this connection, sev- 
eral parameters are to be considered. 

Regarding the material properties, we consider only 
weakly doped or amorphous semiconductors in the regime 
of nearest-neighbor hopping conductance. Hence, we cal- 
culate absorption by close pairs independent of each 
other. To keep the model adequate we have to require 
that the typical inter-center distance within the pair, rc, 
should be much smaller that the typical distance between 
defect centers, f — (47rnd/3)~^/^. Here is the defect 
concentration. The hopping distance rc is discussed in 
Sec. Ifvl. 



Another requirement is that the impurities are not too 
shallow, so that the electrons cannot be excited from the 
localized states to the conduction band by the AC per- 
turbing potential. There are experimental examples of 
this, where suck a excitation serves as a source of non- 
linear behavior.l3 

According to the present calculation, the most interest- 
ing effects occur at "low" frequencies when LOTn{T) <C 1. 
This requirement also ensures that the relaxation absorp- 
tion dominates the resonant one. Certainly, the minimal 
relaxation time t„ (T) is a material property. Usually the 
above inequality is met at low temperatures for frequen- 
cies in the range 100 - 1000 MHz. 

The TLS model in glasses is restricted to very low 
temperatures where higher energy levels are not excited. 
The situation is a bit different in semiconductor materi- 
als where the inter-level splittings are of the order of the 
Bohr energy. Consequently, the nearest-neighbor hop- 
ping conductance can be effective in the temperature 
range up to a few K. 

The main objectives of this paper is to show that 
nonlinear effects are anomalously large. Indeed, in the 
case of deformational absorption an additional parameter 
(tjra)"^/^ ^ 1 is present in the nonlinear expansion (|34|), 
while the the case of piezoelectric interaction the nonlin- 
ear expansion starts from the fist power of dimensionless 
amplitude \F\, Eq, (41). Furthermore, the nonlinear con- 
tributions have pronounced magnetic field dependences 
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different from the linear ones. We hope that those fea- 
tures will allow experimentalists to detect the nonlinear 
behavior and to discriminate between different relaxation 
mechanisms for localized states. 

In course of the present calculations we have assumed 
the inequality e/ Er ^ 1 to he met, see Eq.(12). It is im- 
portant that for e one has to substitute the value which 
gives the dominating contribution in the final integration 
over the pair distribution function. This value is actually 
different for the linear and the nonlinear contributions, 
and it is also dependent on the relaxation mechanism. 
It turns out that for the linear absorption this typical 
e ^ fcsT, while for the nonlinear contributions it is re- 
duced by a factor y/ujTs for the deformational interaction, 
by \F\ for the piezoelectric one, both calculated under 
the condition e ^ Er- A similar estimate is necessary 
to choose a proper value for r in the expression for Er- 
This value depends on the quantities ^furf^, Aq, as well 
as on the magnetic field. Thus the experimental vari- 
ables intensity, frequency and magnetic field, in addition 
to the system parameters tq and Ag, influence the be- 
havior of the nonlinear absorption. This rich parameter 
space allows for a large range of experiments. 



VII. CONCLUSIONS 

In this paper, we have analyzed nonlinear contributions 
to the acoustic and electromagnetic absorption by local- 
ized electron states in semiconductors in the regime of 
hopping conductance. The most important conclusions 
are the following. 

• The total behavior of absorption is determined by 
the product ut where t{T) is the minimal relax- 
ation time for a pair with energy splitting of the 
order fc^T. 

• The anomalous nonlinear behavior occurs at lot <ti 
1. In the case of deformational relaxation mecha- 
nism for the localized electrons a large additional 
factor (wr)"^/^ appears in front of the item oc F^ in 
the expansion of nonlinear absorption. In the case 
of piezoelectric relaxation mechanism the expan- 
sion starts with the item cx |F| rather than cx F^. 
Here F is the dimensional AC field amplitude. 

• The anomalous nonlinear absorption is strongly in- 
fluenced by an external magnetic field, the influence 
being dependent both on the electron pair distribu- 
tion function, on the dominating relaxation mech- 
anism for the localized electrons, and on the direc- 
tion of the field polarization vector with respect to 
the magnetic field. The infiuence of magnetic field 
on the linear absorption and nonlinear corrections 
is substantially different. 

As a result, the physical picture of weakly nonlinear ab- 
sorption appears rich and informative. Our estimates 



show that the effects under consideration are accessible 
for the modern experiment, and many important charac- 
teristics - the dominating relaxation mechanism, the im- 
portance of the Coulomb gap, typical hopping distances, 
etc. - can be extracted by comparison to the present 
theory provided the experiment will be done. 

It should be emphasized that there is a close similar- 
ity between the present and tho-results of our previous 
calculations for glassy materials.u However, the localized 
states in disordered semiconductors, being charged, can 
be influenced by magnetic field which makes them easier 
to investigate. We therefore also hope that the studies of 
semiconductor systems will also provide a new informa- 
tion regarding nonlinear response of TLSs in glasses. 
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